Draft version February 1, 2008 

Preprint typeset using LATp^X style emulateapj v. 16/07/00 



O 
O 

(N 

u 

< 
in 



> 
o 

o 



or 

6 



13 



VLA OBSERVATIONS OF THE H92a LINE FROM NGC 5253 AND HENIZE 2-10 : IONIZED GAS 

AROUND SUPER STAR CLUSTERS 

Niruj R. Mohan 

Raman Research Institute, Bangalore, India 
and 

Indian Institute of Science, Bangalore, India 
AND 

K.R. Anantharamaiah 

Raman Research Institute, Bangalore, India 
AND 

W.M. Goss 

National Radio Astronomy Observatory, Socorro, USA 
Draft version February 1, 2008 

ABSTRACT 

We have detected the H92a radio recombination line from two dwarf starburst galaxies, NGC 5253 
and He 2-10, using the Very Large Array. Both the line data as well as the radio continuum data arc used 
to model the properties of the ionized gas in the centers of these galaxies. We consider a multi-density 
model for radio recombination lines and show why previous models, which were based on the assumption 
of gas at a single density, are valid in many situations. The models show that the ionized gas has a 
density of ^10 4 cm -3 in both galaxies, with an effective size of 2-10 pc and a total mass of ~10 4 Mq. 
The derived production rate of Lyman continuum photons is ~ 2.5 x 10 52 s _1 in both the galaxies and the 
corresponding mass of stars (assuming a Salpeter IMF) is ~ 10 5 Mq. The implied stellar density shows 
that the observed radio recombination lines arise from ionized gas around super star clusters (SSCs) in 
both galaxies (these SSCs have been recently detected through their radio continuum emission). The 
existence of ~10 4 Mq of ionized gas within a few parsecs of an SSC places strict constraints on dynamical 
models. Using simple arguments, the parameter space for a few possible models are derived. The well 
known radio-FIR correlation also holds for NGC 5253, although the radio emission from this galaxy is 
almost completely thermal. It is shown that NGC 5253 is strong evidence that the component of FIR 
emission from warm dust is correlated separately with the component of radio emission from thermal 
bremsstrahlung. 

Subject headings: galaxies: individual (NGC 5253, He 2-10) — galaxies: ISM — galaxies: starburst — 
galaxies: star clusters — radio lines: galaxies 



1. INTRODUCTION 

It is now known that star formation in starburst galax- 
ies is not just a scaled-up version of quiescent star forma- 
tion in normal galaxies. For example, the phenomenon 
of starbursts involves the formation of super star clusters 
(SSCs), which are believed to be the pre-cursors of the 
present-day globular clusters (see recent review by Whit- 
more (2000) and references therein). These SSCs may gen- 
erate large-scale superwind outflows (Heckman, Armus, 
& Miley 1990) which can carry away a good fraction of 
the gas and processed heavy elements from the starburst 
region. Star formation also occurs in the form of indi- 
vidual field stars and small unresolved clusters. Along 
with stellar emission diagnostics, studies of ionized gas 
around OB stars provide information on the young stel- 
lar population and also characterizes the surrounding in- 
terstellar medium. Radio continuum and recombination 
line emissions are good extinction-free probes of this gas. 
High resolution radio continuum imaging of starburst cen- 
ters, for example, has led to the discovery of a host of 
pc-sized thermal and non-thermal compact sources (NGC 
253: Ulvestad & Antonucci (1997) and references therein; 



NGC 2146: Tarchi et al. (2000); M 82: Kronberg, Bier- 
mann, & Schwab (1985); NGC 4038/9: Neff & Ulvestad 
(2000)). Studying the compact thermal sources is difficult 
in the UV, optical or even in IR, owing to extinction prob- 
lems. A nebular diagnostic like radio recombination lines 
(RRLs) can probe thermal gas with large emission mea- 
sures. In the cm- wavelength range, RRLs have been shown 
to be detectable from nuclear starburst regions in a num- 
ber of nearby galaxies (Shaver, Churchwell, & Rots 1977; 
Seaquist & Bell 1977; Puxley et al. 1991; Anantharamaiah 
et al. 1993; Zhao et al. 1996; Phookun, Anantharamaiah, 
& Goss 1998). These lines have proven to be good probes 
of gas density, and are reliable estimators of ionized gas 
mass and the ionizing photon flux in these starburst re- 
gions. RRLs have been studied in about 14 galaxies so far 
(Anantharamaiah et al. (2000) and references therein). In 
this paper we present observations of the H92a RRL at 8.3 
GHz from two prototypical Blue Compact Dwarf Galaxies 
(BCDGs), NGC 5253 and He 2-10. 

BCDGs are dwarf galaxies which are dominated by an 
ongoing starburst, with or without a background quies- 
cent star formation. NGC 5253 and He 2-10 are classic 
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examples (Davies, Sugai, & Ward 1998) of BCDGs. NGC 
5253 (D=4.1 Mpc, Saha et al. 1995) is a dwarf galaxy 
with an underlying elliptical structure (Caldwell & Phillips 
1989). This galaxy has large scale ionized filaments (Gra- 
ham 1981) and a number of young star clusters in the halo 
(van den Bergh 1980). The central kpc region is under- 
going active star formation and hosts a number of SSCs 
(Turner, Beck, & Ho 2000; Calzetti et al. 1997; Gorjian 
1996; Beck et al. 1996; Meurer et al. 1995; Caldwell & 
Phillips 1989). Wolf-Rayet emission has been detected in 
the central region (Campbell, Terlevich, & Melnick 1986; 
Walsh & Roy 1987, 1989; Schaerer et al. 1997). Based on 
modeling of the stellar emission and through the presence 
of ambient ionized gas and Wolf-Rayet features (Gorjian 
1996), a spatial age gradient in the star formation in the 
core has been inferred. NGC 5253 is notable for having a 
flat radio continuum spectrum (Lequeux 1971; Beck et al. 
1996), which is attributed to lack of non-thermal radiation 
due to the absence of supernova remnants. This absence 
is believed to be due to the very young age of the star for- 
mation activity. The resulting age estimate is consistent 
with that derived from optical and IR spectrophotometry 
(Beck et al. 1996). The Ha structure in the inner 20" is 
complex and there is evidence for a patchy dust distribu- 
tion (Calzetti et al. 1997; hereafter CMB97) . Recently, 
Turner, Beck, & Ho (2000), hereafter TBH00, have made 
a high resolution image of the central region in the radio 
continuum at 15 and 22 GHz and have detected a compact 
source which they interpret as a 'supernebula' surrounding 
a nascent SSC. 

He 2-10 (D=9 Mpc) was the first Wolf-Rayet galaxy 
discovered (Allen, Wright & Goss 1976). Three star form- 
ing regions have been identified in the central part of the 
galaxy and among them, the central ~5" region, called re- 
gion 'A', is the youngest (Mendez et al. 1999). This region 
exhibits predominant Wolf-Rayet spectral features and is 
also the most active star forming region. HST imaging in 
the UV has shown the presence of a series of young clus- 
ters, believed to be SSCs (Conti & Vacca 1994; Johnson et 
al. 2000). Deep Ha images of this galaxy (Mendez et al. 
1999) indicate a large scale bipolar outflow of gas. Molec- 
ular gas, imaged through CO lines, shows a disturbed 
morphology with a rotating disk or a bar (Kobulnicky et 
al. 1995). Radio continuum from He 2-10 has an overall 
spectral index of —0.6 and is therefore dominated by non- 
thermal emission (Allen, Wright & Goss 1976). Recent 
high resolution radio continuum imaging by Kobulnicky 
& Johnson (1999), hereafter KJ99, has shown the exis- 
tence of five compact radio knots with flat spectra. These 
radio knots are interpreted by the authors as high den- 
sity nebulae surrounding young SSCs, which are probably 
totally obscured in the optical. 

Based on their nuclear Br7 flux densities and the pres- 
ence of active star formation in the central region, we se- 
lected these two galaxies to search for RRLs from the cen- 
tral few hundred parsecs. Our observations are a part of a 
program we have undertaken to detect RRLs in nearby 
starbursts and to follow up the detected sources using 
multi- frequency RRL and continuum observations. The 
aim is to model the density structure, kinematics and ge- 
ometry of the ionized gas in these regions and to deter- 
mine an extinction-independent star formation rate. An 



example of such a study of Arp 220 has been presented by 
Anantharamaiah et al. (2000). 

2. OBSERVATIONS AND RESULTS 

NGC 5253 and He 2-10 were observed with the Very 
Large Array (VLA) in both the CnB and the DnC con- 
figurations with an on-source integration time of 3 hours 
per source in each configuration. Some details of the ob- 
servations are given in Table 1. The system was tuned to 
detect the H92a line {v rest =8309.385 MHz) using a total 
bandwidth of 25 MHz and 15 spectral channels. Data from 
both polarization channels were recorded. The data was 
hanning-smoothed offline to reduce the effects of Gibbs 
ringing. The effective velocity resolution is 113 km s _1 . 
The flux density scale was fixed by using 3C 286 as a cali- 
brator and phase and bandpass calibrations were done by 
using a nearby point source. All the data analysis were 
done by using standard algorithms within the software 
AIPS. The continuum visibility data was self-calibrated 
and the resultant antenna gains were applied to the line 
data as well. The CnB and DnC array data of each galaxy 
were later concatenated and combined (C+D) continuum 
and line datasets (with a resolution of ~5") were obtained. 
The final line images were derived after subtracting the 
continuum using the AIPS task UVLIN which subtracts a 
linear function of frequency from the real and imaginary 
parts of the observed visibilities (Cornwell, Uson & Had- 
dad 1992). The continuum and line data are presented in 
Tabic 2. Figure 1 shows the continuum images of the two 
galaxies. The images show evidence for filaments, espe- 
cially in the case of NGC 5253, where there are extensions 
in the southern direction. These filaments are also seen 
in optical line and continuum images and extend out to a 
few arcminutes. The total continuum flux density of NGC 
5253 at 8.3 GHz is 58 mJy. The radio spectrum of this 
galaxy is known to be flat from 1.4 GHz up to 23 GHz, 
based on both single-dish and scaled array VLA observa- 
tions (see Lequeux (1971) and also Turner, Ho & Beck 
(1998) and references therein). Lower resolution images 
indicate that there is emission over at least 2'. The short- 
est spacing available in the D configuration of the VLA 
implies that our images are not sensitive to structures on 
scales larger than ~2'. We estimate that the total flux den- 
sity of NGC 5253 in our image is 80% of the single-dish 
flux density at 8.3 GHz. The radio continuum emission in 
NGC 5253 peaks around region 'A' (the cm-emission peak 
in Turner, Ho and Beck 1998) which contains the two main 
star forming regions - clusters 4 and 5 (CMB97) . The ra- 
dio continuum emission is weak to the south of the peak, 
where the older clusters are believed to have dispersed the 
gas. The radio emission in He 2-10 covers the two main 
star forming regions - 'A' and 'B' (region 'A' is the cen- 
tral 6" region and region 'B' is the secondary star forming 
region 8.5" east of 'A', Corbin, Korista & Vacca (1993)). 
Figure 2 shows the H92a line spectra of the two galaxies 
at the position of the peak continuum emission. The line 
emission in both the galaxies is spatially unresolved, even 
in the higher resolution C array data. 

3. MODELING THE H92a LINE 

In this section, we model the observed data in terms of 
the local electron density, temperature and size of the ion- 
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ized gas. The models for the ionized gas are based on the 
observed H92a line flux density and the continuum flux 
densities at 8.3 GHz, 4.8 GHz and 15 GHz (see Table 2) 
from the region of RRL emission. Since the line forma- 
tion process is dependent on the geometry of the emitting 
gas, we consider four different models. Model I is for a 
single spherical HII region of fixed size and constant den- 
sity. Model II consists of a uniform slab of constant den- 
sity, and fixed thickness and beam filling factor. Model III 
consists of a collection of HII regions, of fixed density and 
size. The number of such HII regions is normalized using 
the observed line flux density. The third model is sim- 
ilar to those described by Anantharamaiah et al. (1993) 
and further adapted by Zhao et al. (1996) and Phookun, 
Anantharamaiah, & Goss (1998) for modeling RRLs from 
ionized gas in starburst regions. For all three models, a 
parameter search is done in density n e in the range 1CP 2 
cm~ 3 < n e < 10 6 cm~ 3 and size / in the range 0.01 pc < 
I < L pc, where L is the upper limit to the size of the line 
emitting region given by the beam size, which is 150 pc for 
NGC 5253 and 315 pc for He 2-10. Model II also includes 
the beam filling factor as a free parameter, which is the 
lateral extent of the slab of ionized gas within the region 
of total lateral size L pc. We also consider a fourth model 
(Model IV) in which the co-existence of multiple densi- 
ties is assumed. In this model, we assume a population 
of HII regions with different densities, but with each HII 
region having a constant density. We consider those den- 
sities which correspond to compact (n e > 10 3 cm~ 3 and 
diameter d between 0.1 and 1 pc) and ultracompact HII 
regions (n e > 10 4 cm~ 3 and d < 0.1 pc). A HII region 
of density n e is assigned a size I given by the relation n e 
= Ci l a . The parameters a and Ci depend on the ambient 
medium and the evolutionary phase of the HII region. In 
our model, these are free parameters. There is evidence 
for the existence of such a power law for galactic compact 
and ultracompact HII regions (see, for example, Habing 
and Israel 1979 or Garay et al. 1993). The population of 
HII regions of various densities is described by a power-law 
distribution function, N(n e ) cx nf , which is normalized by 
the observed line flux density. The distribution function 
for galactic diffuse HII regions is a power law (Kennicutt, 
Edgar & Hodge 1989). Comcron & Torra (1996) found 
evidence for a power law luminosity function for galactic 
ultracompact HII regions, although Casassus et al. (2000) 
claim a peaked distribution function. The actual form of 
the distribution would depend in detail on the star for- 
mation history of the starburst region as well as on the 
evolution of the HII regions. In the model, we use (3 as a 
free parameter. The fourth free parameter in Model IV is 
n™", the minimum density in a given distribution. 

The observed data was modeled using models I through 
IV. Valid solutions arc identified as those combinations 
of free parameters which produce both the observed H92a 
line as well as free- free emission at 4.8, 8.3 and 15 GHz less 
than the corresponding total observed continuum flux den- 
sities (which would include both thermal and non-thermal 
components). The gas mass and the ionizing photon flux 
were also derived for each valid solution. All solutions were 
computed for T e = 10000 and 12500 K. These high tem- 
peratures have been derived from optical spectroscopy and 
are consistent with the low metallicity of these two galax- 



ies (Campbell, Terlevich, & Melnick 1986; Walsh & Roy 
1989; Kobulnicky et al. 1997). Model results for the two 
galaxies are summarized in Tables 3 and 4. Since Mod- 
els I and II give similar results, they are grouped together. 
The allowed range of values of densities and sizes of ionized 
gas are listed. The range of densities and sizes allowed by 
the models are quite limited. The gas properties for the 
models with the minimum allowed value of N^ are given 
in the lower half of Tables 3 and 4. The significance of 
the minimum allowed N^c is as follows : For an arbi- 
trary continuum optical depth r c , Nlj, c oc S// x 1 _^1 Te , 
where S// is the free-free radio continuum flux density. 
Hence, for an observed value of S//, the calculated value 
of the required ionizing photon flux Ni yc increases with 
increasing continuum optical depth, t c . A given value of 
Sff thus implies a lower limit to the value of N iyc . We 
take the minimum consistent Nlj, c corresponding to the 
minimum consistent r c for modeling. The values for the 
derived parameters given in Tables 3 and 4 correspond to 
the average of the range of values allowed by the n e -l so- 
lutions. There are two classes of models allowed by the 
existing data : n e < 2500 cm~ 3 , I ^5 pc and n e > 2500 
cm~ 3 and I ^30 pc The model solutions are indicated in 
separate rows for each class of solutions. The predictions 
of some of the valid models for the free- free continuum and 
line flux densities as a function of frequency are shown in 
Figure 3. The calculated values of the free-free continuum 
flux densities are less than the observed total continuum 
flux densities (which are marked as stars in the figure). 
The remaining continuum emission which is unaccounted 
for in our models will have contributions both from ther- 
mal continuum from any non-RRL emitting gas and also 
from the non-thermal continuum (the latter, only in the 
case of He 2-10, since NGC 5253 has very little contribu- 
tion from non-thermal emission in the central region). If 
a significant amount of such a non-RRL emitting thermal 
gas is present, the amount of such gas will be strongly con- 
strained by the continuum flux densities at v < 1 GHz. 
All the above models will be referred to as 'low-resolution' 
models since these models are based on data with a res- 
olution of 5", corresponding to a linear size of 150 pc for 
NGC 5253 and 315 pc for He 2-10. In all four models, the 
favored parameters are densities of n e ~10 4 cm -3 , effec- 
tive sizes of 3-8 pc and masses of ionized gas of ~10 4 M Q 

Tables 3 and 4 show that for both the galaxies, the 
four models give results which are consistent. But the 
minimum value of N^c varies slightly among the various 
models. This variation is due to the different geometries 
assumed in the models and illustrates the uncertainitics 
associated with modeling data with a limited angular res- 
olution. In section 4, we discuss the implications of these 
results of the low-resolution models. 

3.1. Model IV and the single- density approximation 

As is well known, radio recombination line formation is a 
non-LTE process. This property can be used to derive den- 
sities from observations if the geometry is known. From 
the theory of RRLs, it is known that a particular recombi- 
nation line is most efficiently emitted from gas at a particu- 
lar density and that gas at lower densities emits RRLs pre- 
ferentially at lower frequencies and vice versa (Zhao et al. 
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1996). This property of RRLs (as density filters) implies 
that a particular recombination line would receive max- 
imum contributions from a restricted range of gas densi- 
ties. Hence previous work which modeled single RRL data 
using the single-density approximation (corresponding to 
Models I, II and III in this work) were reasonably success- 
ful (Anantharamaiah et al. 1993; Zhao et al. 1996, 1997; 
Phookun, Anantharamaiah, & Goss 1998; Anantharama- 
iah et al. 2000). Therefore, with multi- frequency RRL 
data, the gas can be modeled as a combination of a mini- 
mum number of distinct density components. This method 
was successfully employed to derive a three-density com- 
ponent model of the ionized gas in Arp 220 based on line 
and continuum data from 1.4 GHz to 230 GHz (Anan- 
tharamaiah et al. 2000). 

However, the single-density approximation based on the 
above argument is valid only for certain functional depen- 
dences of emission measure on density (i.e., on the relative 
values of emission measures for gas at different densities : 
the line flux density is most strongly peaked as a function 
of density when f3+2+3/a—0. For the line intensity, the 
condition is /3+2+l/a=0. See section 3 for definitions of 
a and (3 ). For example, Figure 4(a) shows the depen- 
dence of the expected line intensity at 8.3 GHz as a func- 
tion of density, for a constant EM and the line strength 
does indeed have a peak as a function of electron density. 
On the other hand, if a statistically well-sampled popula- 
tion of HII regions with a power law distribution of EM 
with density (as in model IV) is assumed, the line strength 
no longer peaks for a particular density, but instead, is a 
monotonic function of density. For example, it can be 
seen in Figure 4(b) that when N^ yc is held constant (i.c, 

EM oc n^ 3 ), the line intensity is no longer a peaked 
function of density. The results of model IV show that so- 
lutions exist only for values of n™ m ranging between 2500 
and 25000 cm -3 . Also, the solutions with the minimum 
values of Nl vc (the minimum value of Nl vc which is con- 
sistent with the given data, corresponding to the minimum 
needed continuum optical depth, see section 3 for an ex- 
planation) correspond to densities between 5000 and 10000 
cm' 3 . Though no upper limit to the density is imposed 
on the models, the nature of the allowed solutions is such 
that only gas of densities very close to n™" contributes 
to the line flux density and also absorbs almost all of the 
required ionizing UV photons, i.e., all the higher density 
gas in the distribution contributes negligible amounts of 
RRL emission; this result is common to all valid solutions 
in this model. We find that the allowed combinations of 
the free parameters a and \ m i n (the latter being derived 
from n™" and Ci : see section 3) are such that the value of 
n™" 1 is about 10 4 cm' 3 and the size of the gas at this den- 
sity (which is the dominant contributor to the total line 
strength) is between 3-10 pc. This result explains why 
the results of models I, II and III agree well with those 
of model IV. Hence model IV, which is the first attempt 
at a consistent multi-density model for extragalactic RRL 
data, validates the assumption of a single-density ionized 
gas in modeling the line emission. It should be noted that 
the narrow range of allowed densities near ~10 4 cm' 3 is 
not the density at which the H92a line emission is ex- 
pected to peak (from Figure 4(a), this density is ~10 2 
cm -3 ). That both the constant density models as well as 



the multiple density models give the same density of ~10 
cm' 3 , vastly different from the density at which H92a is 
expected to peak, is reason enough to believe that most of 
the ionized gas in this region does indeed have a density 
close to 10 4 cm -3 and that the amount of gas at a density 
of 10 2 cm' 3 must be insignificant compared to the higher 
density gas. 

However, the above argument works only if the density 
distribution is a power law. In general in a multi-density 
environment, where the total EM of gas of a given density 
is an arbitrary function of density, an approach based on 
Model IV will not work. In such a case, multi-frequency 
RRL data is essential and models have to be constructed 
to explain the data with a minimum number of different 
density components (as was done in the case of Arp 220 
by Anantharamaiah et al. 2000). 

4. RESULTS FROM THE LOW- RESOLUTION RRL DATA 

4.1. NGC 5253 

Most of the models discussed in the previous section 
yield electron densities in excess of a few 1000 cm -3 within 
an effective size of 2-10 pc (Tables 3 and 4). The min- 
imum value of the derived Nlj, c is 2.5 xlO 52 s _1 . This 
value of Ni yc is computed from the low-resolution data 
and is consistent with the ionizing photon flux computed 
from other nebular diagnostics listed below. The contin- 
uum flux density of 28 mJy measured from the line emit- 
ting region (which is the position of the peak continuum 
emission, corresponding to the 6.5" x3. 3" sized resolution 
element) in the 8.3 GHz continuum image (Table 2) cor- 
responds to Nl V c — 3xl0 52 s _1 , based on the assumption 
of an optically thin gas. The Br7 flux from the central 5" 
region (Davies, Sugai, & Ward 1998) implies an Nl vc of 
1.5xl0 52 s _1 . The Ha flux density inside an aperture of 
5" (measured using the Ha image of CMB97 ) gives an 
N Lyc of l.OxlO 52 s" 1 . 

Taking the average value of Nl vc = 2.5xl0 52 s' 1 and 
using a Salpeter IMF with m upper =80 Mq and the tables 
given in Vacca, Garmany, & Shull (1996), the number of O 
stars (03-09) is estimated to be ~1800. It should be noted 
that we have used the numbers from Vacca, Garmany, & 
Shull (1996) for solar metallicity whereas the two galaxies 
NGC 5253 and He 2-10 have sub-solar metallicities; there- 
fore the calculated values of Ni yc will be underestimated. 
However, this error is less than 20% (Vacca 1994). As- 
suming mi ower to be 1.0 M©, the total number of stars 
within a region of size ^5 pc is ^3.5xl0 5 . This ^5 pc 
sized region of gas could be spread over an entire resolu- 
tion element (150 pc for NGC 5253). On the other hand, 
if the gas is physically in a compact region of size ~5 pc, 
then the corresponding stellar density indicates that the 
ionizing sources may well be SSCs (Meurer et al. 1995). 

4.2. Henize 2-10 

Similar to the results for NGC 5253, most models for 
He 2-10 yield electron densities in excess of a few 1000 
cm' 3 within an effective size of 2-10 pc (Table 4). Gas 
densities of 1000 cm' 3 have been derived in the central 
2" region using [SII] line ratios (Sugai & Taniguchi 1992). 
The minimum derived value of Ni yc for the 5" sized RRL 
emitting region is 4xl0 52 s _1 . The central 2" region has 
been resolved into a number of star-forming knots and 
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the value of ~NLyc derived from observations which resolve 
these knots is comparable to that derived here. Johnson et 
al. (2000) obtain a value of N Lyc = 2.5xl0 52 s" 1 from their 
UV continuum images, while Doyon, Puxley, & Joseph 
(1992) derive a value of 4xl0 52 s _1 from their measured 
Br7 flux inside a 5" aperture. For an Nl V c of 4xl0 52 s -1 , 
the number of O stars required is ^2900 and the total 
stellar mass is ^5.5xl0 5 M (using a Salpeter IMF with 
m up per=80 Mq and mi ower to be 1.0 M , and the results 
of Vacca, Garmany, & Shull 1996). As in the case of NGC 
5253, the corresponding stellar density for a 3-10 pc region 
corresponds to an SSC, if we assume that the ionized gas 
is physically confined to this ~5 pc region (see section 5.2 
for more details). 

Though the derived ionizing photon flux and the length- 
scale of the gas for both galaxies imply an SSC, the 8.3 
GHz RRL data alone does not clarify whether the ion- 
ized gas is in a single clump of diameter ~5 pc or it is 
spread out in some unknown geometry over the entire 5" 
region. In the latter case, the ionizing source need not be 
an SSC. In the next section, we present arguments that the 
ionized gas is indeed much less than 5" in size. We also 
discuss the implications of combining our RRL data and 
model results with recently published sub-arcsecond reso- 
lution continuum images of the centers of the two galaxies 
(obtained by TBH00 for NGC 5253 and by KJ99 for He 
2-10). 

5. COMPACTNESS OF THE RRL EMISSION : EVIDENCE 
FROM OTHER HIGH RESOLUTION OBSERVATIONS 

5.1. NGC 5253 

We can indirectly infer that the RRL emitting region is 
smaller than 5" from the Ha and H/3 images of CMB97 . 
The Ha image in the inner 5" is quite clumpy with exten- 
sive filaments and knots. Though the RRL emission need 
not follow the Ha emission, it is unlikely that the high den- 
sity gas responsible for RRL emission is spread smoothly 
over the 150 pc region, given the patchy appearance of the 
Ha emission, and hence we may conclude that the RRL 
emitting gas must also be patchy within the 5" beam. 

TBH00 have imaged NGC 5253 at 23 GHz and 15 GHz 
with the VLA with angular resolutions of 0.13" and 0.2", 
respectively. A compact (~0.1", or 2 pc) source with a 
flux density of 11 mJy at 15 GHz, and hence a brightness 
temperature of 12000 ± 3000 K, was discovered. This was 
interpreted as a HII region, and in addition, these authors 
detect radio emission at a low level with a size of 1-2". 
Further confirmation of the thermal nature of this source 
comes from recently detected IR emission from the cen- 
tral object at 11.7 and 18.7 ^m (Gorjian, Turner, & Beck 
2001). Though the RRL emission is not expected to ex- 
actly follow the morphology of continuum emission (since 
a particular RRL comes from a limited range of densi- 
ties, and the continuum emission is sensitive to all densi- 
ties if optically thin), no line emission is expected from 
continuum-free regions and hence we can associate the 
observed H92a emission with the inner 1-2" region. As- 
suming that the observed radio continuum emission from 
this galaxy is purely due to free-free emission, TBH00 
derive an rms electron density of 4 x 10 4 cm~ 3 and an 
N_Ly C of 4 x 10 52 s _1 . From the ionization requirement 
and the size of the object, TBH00 suggest that this com- 



pact source is a 'supcrnebula' surrounding an SSC. Since 
this source is the strongest radio continuum source in the 
inner 5" region and since the derived properties of the 
'supernebula', especially the local electron density, match 
those derived from the RRL data, we tentatively identify 
the RRL emitting gas with the continuum source detected 
by TBH00 . Therefore, we have modeled the 'superneb- 
ula' to estimate the fraction of the observed RRL emission 
this compact source is capable of producing. The models 
are constructed to produce the observed 11 mJy contin- 
uum emission at 15 GHz from gas of size ~0.1". In addi- 
tion, the continuum and RRL emission by this gas at 8.3 
GHz is also constrained to be less than the observed val- 
ues (the latter two constraints are imposed as upper lim- 
its since the observed values correspond to the larger 5" 
area). The results of this model are summarized in Table 
5. Only 10-35% of the observed 8.3 GHz line emission can 
be explained by the TBH00 'supernebula'. All densities > 
5xl0 4 cm -3 are consistent with the existing observations. 
The remaining line emission presumably originates in the 
surrounding larger nebulosity with a size of 1-2" . A calcu- 
lation similar to the one above, where we model the inner 
3" region with a 15 GHz flux density of 22 mJy (TBH00), 
shows that such a model is indeed capable of explaining 
the entire observed RRL line emission. 

The 'supernebula', as imaged by TBH00 is separated by 
about 0.9" from the position of the brightest Ha knot in 
CMB97 . This knot was identified with the UV Cluster 5 in 
the HST image by Meurer et al. (1995). Since the la error 
between the absolute reference frames of the Guide Star 
Catalogue 1 of the HST, on which the Ha and UV images 
are based, and the extragalactic radio reference frame, on 
which the images in TBH00 are based, is 0.5" (Russell et 
al 1990 and also Calzetti, and Rajagopal, private commu- 
nication), the radio 'supernebula' of TBH00 can also be 
identified with the Ha peak (and hence the UV cluster 5) 
within these errors. The significance of this identification 
is two-fold. (1) An additional highly obscured SSC need 
not be invoked and (2) we can relate the known properties 
of the optical cluster with those of the 'supernebula'. The 
Ha flux for this cluster is 4.8 xl0~ 16 W m~ 2 (sec CMB97 ) 
from which we infer the value of Ni yc to be ~1.3x 10 51 s _1 . 
The discrepancy between this value and N^ >2.5xl0 52 
s _1 which is derived from radio continuum data (Table 
5) can be attributed to extinction. However, as the age of 
the optical cluster is well constrained, the dynamical age of 
the ionized 'supernebula' around it can also be constrained 
(see section 7.1 for further details). 

5.2. Henize 2-10 

The central 5" region harbors about nine UV clusters, 
possibly SSCs, arranged in an arc of size ~2" (Conti & 
Vacca 1994). This region has a large patchy extinction 
and also clumpy ionized gas. The Br7 , UV and optical 
continuum images of the center show the presence of sev- 
eral knots (Davies, Sugai, & Ward 1998; Conti & Vacca 
1994; Johnson et al. 2000). KJ99 have imaged the radio 
continuum emission of the nucleus of Henize 2-10 with an 
angular resolution of 0.5" using the VLA. They detect 4-5 
compact knots, each of diameter ~0.6", arranged within 
a 5"x2" region. These knots of star formation detected 
individually in the optical, Ha , Br7 and radio continuum 
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are not positionally consistent, possibly due to intrinsic as- 
trometric errors. Attempts at aligning the various images 
using global astrometric offset corrections have shown that 
not all knots have counterparts at all wavelengths (KJ99 
and Johnson et al. 2000); thus there appears to be hidden 
regions of star formation in the central 5" region. This 
morphological evidence for the clumpy gas at the center 
indicates that the RRL emitting gas is also clumpy on 
scales <5" . 

Based on the continuum spectra, KJ99 infer that the 
radio knots consist of thermal gas at densities 1500-5000 
cm~ 3 (optically thick at 8.3 GHz) and also derive a total 
N_Ly C of 4x 10 52 s _1 . They suggest that each compact knot 
may be ionized by a SSC. We have modeled this region to 
look for solutions which produce the observed radio contin- 
uum flux densities at 15, 8 and 5 GHz, as reported by KJ99 
. The allowed solutions are also constrained to produce ei- 
ther all or a part of the observed H92a line. The model 
results are listed in Table 5 and these radio knots, if ther- 
mal, are indeed capable of producing all of the observed 
H92a recombination line. Therefore we identify these mul- 
tiple radio continuum knots with the gas responsible for 
emitting the observed RRL at 8.3 GHz and associate the 
ionizing sources with SSCs. 

6. NGC 5253 AND RADIO-FIR CORRELATION 

The FIR luminosity and the radio luminosity of galax- 
ies which are dominated by star formation as opposed to 
an AGN are very tightly correlated with an rms of <0.2 
dex on the logarithmic scale (Condon 1992). This radio- 
FIR correlation holds for several types of galaxies such 
as normal spirals, dwarf starbursts and dwarf irregulars 
(see Condon (1992)). The physical reason for this correla- 
tion is not well understood in detail but may result from 
a chain of events involving star formation. The thermal 
(free- free) portion of the radio emission is from gas ionized 
by OB stars and the non-thermal (synchrotron) portion of 
the radio emission arises from the cosmic rays accelerated 
in SNRs, the end point of Type II SNe caused by massive 
stars. The warm (T > 30K, Seo^m/Sioo^m > 1) dust por- 
tion of the FIR emission arises from dust directly heated 
by OB stars and the cool (T < 30 K, S 60tJ .m/Swa^m < 0.8) 
dust portion of the FIR emission arises from dust heated 
by the diffuse interstellar radiation field. Xu et al. (1992) 
have shown in the case of the LMC that the warm dust 
emission correlates with the thermal radio emission, while 
the cool dust emission correlates independently with the 
non-thermal radio emission. Also, the ratio of the two cor- 
related quantities for both these independent correlations 
was shown to be the same as that between the the total 
FIR and the total radio emission strengths. 

The radio continuum of most disk galaxies is domi- 
nated by non-thermal emission since the thermal frac- 
tion is typically 0.05-0.2 at 1.4 GHz (Condon 1992). But 
it is known that the BCDGs have a higher fraction of 
thermal emission than disk galaxies: Klein & Wunder- 
lich (1987) showed that the average radio spectral index 
(a) is —0.4 for BCDGs, while it is —0.7 for normal spi- 
rals, oc v a . Though the mean spectral index is indeed 
flatter for BCDGs, the dispersion in the index is much 
larger than for spirals and hence they did not estimate 
an average thermal fraction for the BCD galaxy sample. 



BCDGs also have a higher dust temperature (Wunder- 
lich and Klein 1988, hereafter WK88) compared to spiral 
galaxies which implies that their FIR emission is domi- 
nated by emission from warm dust. It has been suggested 
by Klein & Wunderlich (1987) and WK88 that the reason 
the BCDGs also obey the radio-FIR correlation is because 
this dominant warm dust emission correlates with the ther- 
mal radio emission, which is the dominant portion of the 
radio continuum. 

NGC 5253 may well be an unique case which will help to 
understand the details of the radio-FIR correlation. Radio 
continuum flux density measurements at various frequen- 
cies using single dishes and low resolution images have 
shown that the radio emission of NGC 5253 is consistent 
with being almost purely thermal (see Beck et al. (1996) 
and also the spectral index map of Turner, Ho, & Beck 
1998). Its total flux density at 1.4 GHz is 89 mJy (from 
the NVSS catalogue of Condon et al. 1998 and single-dish 
observations) and the FIR flux is calculated to be 8xl0 8 
L & from the IRAS flux densities. From the ratio of the 
radio to the FIR luminosity, it can be seen that this galaxy 
obeys the radio-FIR correlation (Condon 1992) to within 
1.5o\ This is a surprising result, since the non-thermal 
emission in this galaxy is less than 10% of the total emis- 
sion at 5 GHz (Beck et al. 1996; Turner, Ho & Beck 1998). 
It is known that the radio flux density and the FIR flux 
density of galactic HII regions are also correlated, but with 
a different ratio (Wynn- Williams & Becklin 1974; Fuerst, 
Reich & Sofue 1987; Broadbent, Osborne & Haslam 1989). 
On the other hand, the ratio of radio to FIR flux densities 
of compact HII regions is not a constant, but depends on 
a variety of factors (Kurtz, Churchwell & Wood 1994). It 
is worth noting that though the radio emission of NGC 
5253 is almost purely thermal, the radio to FIR emission 
ratio is similar to star forming galaxies rather than HII 
regions. One way to explain the correlation is the follow- 
ing : the dust temperature in the nucleus of NGC 5253 is 
high; Calzetti et al. (1995) estimated it to be 51 K. This 
high dust temperature is the reason why the ratio of the 
strength of the warm dust emission to the cool dust emis- 
sion is about 16, one of the largest among BCDGs and disk 
galaxies (Calzetti et al. 1995). Therefore we suggest that 
the reason NGC 5253 obeys the radio-FIR correlation is 
that the warm dust emission and the thermal radio emis- 
sion are indeed correlated with each other, and the slope 
of the correlation is the same as that between the total 
FIR and total radio in normal spirals, as suggested by Xu 
et al. (1992) and WK88 . WK88 had also put forward an 
alternative explanation for BCDGs obeying the radio-FIR 
correlation. They had suggested that a possible lower dust 
fraction in these galaxies would compensate for the lower 
non-thermal radio fraction. This explanation is unlikely 
for NGC 5253, since among all BCDGs, it has the low- 
est non-thermal fraction but the highest dust to gas ratio 
(Lisenfeld & Ferrara 1998). 

7. DISCUSSION 

The detection of a RRL from the central region of 
both galaxies with parameters similar to those derived by 
TBH00 and KJ99 is added confirmation for the thermal 
nature of the central compact source (s). Additionally, an 
RRL study of dense ionized gas has the following advan- 
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tages : (1) Instead of the rms electron density as given 
by the continuum data, local electron densities, indepen- 
dent of filling factor, can be derived from RRL data. (2) 
Therefore, if the line emitting region is spatially resolved, 
then the actual value of Nl vc can be derived (and not just 
the lower limit). (3) Due to non-LTE effects, RRLs retain 
their sensitivity to density even in gas which is moderately 
optically thick, unlike continuum emission. 

7.1. Dynamics of ionized gas around SSCs 

A question not yet addressed regards the dynamics and 
time-scales involved in explaining the existence of ~10 4 
Mq of ionized gas within a few parsecs of an SSC. In this 
section, we will identify the 'supernebula(e)' with the RRL 
emitting gas in both the galaxies and also tentatively asso- 
ciate the 'supernebula' in NGC 5253 with the star cluster 
in the center. The observed FWHM of ~100 km s -1 (es- 
timated from the spectrum before offline hanning smooth- 
ing) within a 2-8 pc sized region implies a dynamical ex- 
pansion time of 2-8 x 10 4 years, quite short compared to the 
estimated age of 1—3 xlO 6 years for the central cluster in 
NGC 5253. Therefore, either the supernebula must have 
been set into motion extremely recently, or there must ex- 
ist a confining mechanism operating in the central region. 
In addition, strong stellar winds from OB and Wolf-Rayet 
stars as well as the natural expansion of the HII region 
due to overpressure would drive the ionized gas outwards. 
Given the large mass, high density and small size of the 
observed HII region (~2 pc for the 'supernebula' of TBH00 
and KJ99 and 5-8 pc for the RRL emitting gas), dynam- 
ical models should constrain both the age of the SSC as 
well as the density of the ambient unshocked medium. For 
NGC 5253, for example, if the 'supernebula' is identified 
with the UV Cluster 5, then the age is determined to be 
less than 3 million years (CMB97). 

If we assume the HII region to have reached a diameter 
of 1 pc solely due to expansion driven by overpressure in 
~2 Myr, then the radius of the initial Stromgren sphere 
would be ^0.01 pc (Spitzer 1968). If at the time of forma- 
tion of the initial Stromgren sphere, the ionizing photon 
flux were ~T0 52 s _1 , then the corresponding initial density 
should be ~2xl0 7 cm~ 3 . In such a case, the mass of the 
shocked neutral gas between the shock and the ionization 
fronts would be more than 100 times greater in the ionized 
gas inside the ionization front (Savedoff & Greene 1955). 

In addition, stellar winds will definitely play an im- 
portant dynamical role. Using StarBurst 99 (Leitherer 
et al. 1999), a mechanical luminosity of 10 38 ergs s _1 is 
estimated for an N^ yc of 10 52 s _1 , assuming a continu- 
ous starburst model. Assuming that the isothermal phase 
sets in early in the lifetime of the HII region, the radius 
of the thin swept-up ionized shell is given by Ki SO = 16 
(^y) 1 ' 4 t,!/ 2 pc, where L is the mechanical luminosity 

in units of 10 36 ergs s _1 , V is the velocity of the stellar 
wind at infinity (assumed to be 1000 km s _1 ), n Q is the 
density of the ambient medium and t$ is the time in mil- 
lions of years (McCray 1983). The thickness of the ionized 
shell is given by n Q Ri SO /(3n s ) (Dyson & Williams 1997), 
where n s , the density of the shocked ionized gas, is given 
by n s =n xV 2 /C 2 for an isothermal shock (V s is the ve- 
locity of the shock, dR iso /dt, and C is the speed of sound 
in the ambient medium) . The values of n Q and t can be ex- 



pressed in terms of n s and the mass of the ionized gas. The 
range of input parameters are 10 3 cm~ 3 < n s < 5xl0 4 
cm" 3 and 10 3 M Q < M < 10 4 M , which are similar to 
the ranges for the models in Tables 3 and 4. The ranges of 
the derived values are 0.01 cm~ 3 < n Q < 100 cm~ 3 and 10 5 
yr < t < 2xl0 6 yr. Although the derived values do not 
violate any known constraints, the allowed solution space 
and a detailed modeling based on future higher resolution 
continuum images would be useful to understand the dy- 
namics and the ambient environment of these regions. For 
an initial attempt at such an analysis, see Tan & McKee 
(2000). 

7.2. High density ionized gas in starburst regions 

The commonly occuring relatively high density of ^10 4 
cm^ 3 for the compact ionized gas in the starburst regions 
in He 2-10 and NGC 5253 and also in other starburst galax- 
ies (Gilbert et al 2000, Zhao ct al 1997, Anantharamaiah et 
al 2000, TBH00, KJ99 ) is striking. One model which may 
be able to explain this commonly occuring high densities 
is that of Jog & Das (1992, 1996). In their model, a GMC 
in the center of a starburst region is shocked by the over- 
pressure of the surrounding molecular inter-cloud medium 
(ICM). The outer shell of the GMC is shock-compressed 
and moves inwards, and at a later stage, due to onset of 
gravitational instabilities, it fragments into star forming 
centers. This model can well explain the formation of SSCs 
of mass >10 5 Mq. An estimate of the density of the shell 
fragments when star formation begins in the imploding 
shell is possible. From Jog & Das (1996), the velocity of 
the shock, V s , at the onset of gravitational instabilities can 
be as high as ~30 km s _1 and the radius of the shock front 
reaches ~40% of the GMC radius. Adopting a value of the 
sound speed, C s , inside the GMC as 4 km s _1 , the den- 
sity of the shocked gas, n s =n xV 2 /C 2 , is increased and 
the value of the density enhancement n s /n can be greater 
than 55 (this enhancement can be up to three times larger 
than this value, for the observed range of values of GMC 
and ICM properties). In their model, Jog & Das (1992, 
1996) assume that the density of the GMC falls radially 
as r _1 and that the density at the boundary is at least 
100 cm~ 3 . Therefore this density enhancement could ex- 
plain the observed densities of 10 3 -10 4 cm~ 3 observed in 
RRLs. It remains to be seen if there is indeed a high pres- 
sure reservoir of molecular gas or shocked hot gas at the 
centers of these galaxies, and if so, whether star formation 
occurs via this process. It should be noted that CO images 
of the centers of these two galaxies do not show evidence 
for gas in this region and therefore we appeal to whichever 
gas reservoir gave rise to the SSCs in the center to provide 
the neccesary high pressure. However the high density am- 
bient medium is not needed as in the previous cases, since 
remnant clouds of density ~10 4 cm~ 3 would still remain 
in the vicinity of the star cluster and be ionized. 

The RRL observations in this paper and the recent con- 
tinuum observations of 'supernebulac' may well be the first 
detection of ionized gas around extremely young SSCs. 
Since these clusters are usually too obscured to be stud- 
ied even in the IR, higher resolution radio continuum and 
line data can be used to infer the properties of the ambi- 
ent medium around SSCs, and ultimately, the formation 
mechanism of these massive clusters. Sensitive line data 
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will also yield kinematic information about these nebulae. 
But an important problem is the identification of the radio 
components with the optical and UV components at the 
sub-arcsecond level, which should be achievable in a few 
years. 

8. SUMMARY 

We have detected the H92a recombination line from two 
dwarf starbursts NGC 5253 and He 2-10. The ionized gas 
responsible for both line emission and the radio continuum 
emission is modeled in terms of electron density and size. 
The line emitting gas in both galaxies is found to be of 
densities ^10 4 cm~ 3 with a size of 2-10 pc, implying an 
Ni yc of a few times 10 52 s _1 . We also consider a multi- 
density model for the RRL-emitting gas and show that 
the assumption of a constant density in previous models is 
valid in a wide range of cases. Though the line emission is 
spatially unresolved and the linear size of the beam area 
corresponds to 150 pc and 315 pc respectively for NGC 
5253 and He 2-10, we present arguments which suggest 
that the line emission is compact (few parsecs). Based on 
the derived parameters, these lines seem to originate from 
ionized gas associated with SSCs in both galaxies. The 
possible 'supernebulae' around these clusters have been 
detected in radio continuum emission and we show that 
these nebulae are capable of producing all the observed 
line emission in both the galaxies. 

Both the detection of RRL emission as well as direct 
radio continuum imaging of these 'supernebula' around 
SSCs pose problems about the lifetimes and dynamics of 
the ionized gas. We present simple arguments to show that 
reconciling the dynamical age of the 'supernebula' with the 



age of the SSCs might be difficult and that detailed obser- 
vational data as well as theoretical modeling are neccesary 
in order to understand the dynamics of HII regions around 
these massive clusters. Finally, we point out that though 
NGC 5253 has a thermal radio spectrum, it seems to 
obey the radio-FIR correlation, which is known to hold for 
galaxies dominated by non-thermal radio emission. This 
feature may well be evidence for the warm dust portion of 
the FIR emission to be correlated with the thermal radio 
emission, independent of the cold dust FIR emission and 
the non-thermal radio emission, as suggested by Xu et al. 
(1992). 
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Fig. 1. — 8.3 GHz radio continuum images of NGC 5253 (left) and He 2-10 (right) observed using the combined C and D array VLA data, 
made with 'natural weighting' of the visibilities. The contours start at 5<x (150/iJy and 160/iJy for NGC 5253 and He 2-10 respectively) and 
increase in steps of V2. The beam sizes arc 6.5"x3.3", PA=-40° for NGC 5253 and 5.0"x4.3", PA=24° for He 2-10. For further details, 
see Table 2. 
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Fig. 2. — The harming smoothed H92o spectra at the continuum peaks of NGC 5253 (left) and He 2-10 (right) at 8.3 GHz for the combined 
C and D array VLA data. The line emission is unresolved in both galaxies. The spatial resolution is 6.5"x3.3" for NGC 5253 and 5.0"x4.3" 
for He 2-10. The velocity resolution of the spectrum is 113 km s -1 . 



Table 1 
Log of VLA observations 



Galaxy 


Array 


Date of 


uv range 


Phase 


Bandpass 








Observation 


kA 


cal 


cal 


km s _1 


NGC 5253 


DnC 


19 Feb 1999 


0.7-42 


1334-127 


1334-127 


404 


NGC 5253 


CnB 


10 Nov 1998 


1.3-150 


1334-127 


1334-127 


404 


Henize 2-10 


DnC 


19 Feb 1999 


0.7-46 


0919-260 


0919-260 


873 


Henize 2-10 


CnB 


10 Nov 1998 


2.0-130 


0919-260 


0919-260 


873 



RRLs from NGC 5253 and He 2-10: Ionized gas around SSCs 



11 



E 



E 

C 

• 

■+-> 
C 

o 
o 

"o 
E 
o 




10 100 




Frequency (GHz) 



Frequency (GHz) 



E 



E 

C 

• — 

c 
o 
o 

"o 

E 
k_ 




~ 3 

E 



X 

3 



C 



O : 



O 
O 



O I 
<L> O 
Q_ 



i 1 1 1 mi| — i 1 1 ii 



TTI| 1 I I I llll| P] 



10 100 



o t_ 
"0.1 




t 

// 



LlillJ 1 



I <3 



10 100 



Frequency (GHz) 



Frequency (GHz) 



Fig. 3. — Predictions for the thermal continuum and the peak RRL flux densities as a function of frequency, for NGC 5253 (top) and He 
2-10 (bottom) for a few valid solutions corresponding to Model I (see section 3 for details). The data points are marked as stars. The solid 
line corresponds to n e =10 3 cm~ 3 , 1=24 pc, T e =10000 K, the dashed line to n e =5000 cm -3 , 1=6 pc, T e =10000 K and the dot-dashed line to 
n e =10 4 cm" 3 , 1=4 pc, T e =10000 K. The predictions of Models II-IV are similar. 
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Fig. 4. — Expected H92« intensity as a function of density for (a) constant EM=10 7 pc cm -6 and (b) constant Nij /C =3xl0 43 s~ x . Pressure 
broadening is neglected and non-LTE effects are taken into account. The curve for constant J^Lyc illustrates the point that for an arbitrary 
population of multi-density ionized gas, the RRL strength as a function of density need not a have a maximum over a restricted range of 
densities. 
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Table 2 
Observational data 



Parameter 


NGC 5253 


He 2-10 


Beam size 


6.5"x3.3" 


5.0"x4.3" 


P.A. of the beam 


-40° 


24° 


Peak continuum flux 






density at 8.3 GHz a 


28 mjy 


13 mjy 


Total continuum flux 






density at 8.3 GHz b 


58 mjy 


29 mjy 


Noise in continuum image (la) 
Peak H92a line flux density 


27 ^Jy 


21 ^Jy 


0.53 mjy 


0.16 mjy 


Line width (FWHM) 


95 ± 10 km s" 1 


112 ± 20 km s _1 


(before hanning smoothing) 






Noise in channel image (la) 


65 fiJy 


45 ^Jy 


4.8 GHz flux density 


30 mJy d 


18 mJy c 


15 GHz flux density 


26 mjy d 


9 mjy c 



a Peak continuum flux density in the central 5" region 

b Total continuum flux density measured over the entire C+D continuum 
image of 30-45" extent : 80% of the single dish flux density (see section 
2) 

Tine flux density in the central beam area of 5" 

d Continuum flux density for the central 5" : Turner, Ho and Beck 1998 

"Continuum flux density for the central 5" , assuming S„ cx v~ os , see 
Kobulnicky and Johnson 1999 for a summary of the spectral index of He 
2-10 
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Table 3 

The low-resolution model results for NGC 5253 



Parameter 


Models I&II 


Model III 


Model IV 


n e (cm- 3 ) a 


500-1000, 


5000- 10 4 


2500-25000 


5000- 10 4 






size (pc) a ' 6 


20-30, 


5-10 


3-7 


2-7 






EM (pc cm- 6 ) a 


2xl0 7 


10 7 -10 8 


~10 8 




2 x10 s 






Mhii (Mq) 


10 5 ,10 4 


10 4 


10 4 



Solutions for minimum Nlj, c 



^Tyc (S^ ) 


1.5xl0 52 


2.5xl0 52 


3.5xl0 52 


n e (cm 3 ) 


5000- 10 4 


5000- 10 4 


5000- 10 4 


size (pc) b 


2-7 


4-7 


3-7 


EM (pc cm" 6 ) 


2 x10 s 


~10 8 


~10 8 


o /Qcont 
>~>th/ ^>obs 


0.3-0.5 


~0.7 


0.3-0.5 


M HI l (Mq) 


10 4 


10 4 


5000- 10 4 


No. O stars c 


1100 


1825 


2600 


Mass of stars (M Q ) d 


2xl0 5 


3.5 xlO 5 


5xl0 5 



a The two rows correspond to two classes of models, both of 
which are consistent with existing observations. The models pre- 
dict different RRL and continuum flux densities at other frequen- 
cies (see Figure 3) 

b 'Effective' sizes are quoted for Models III and IV, i.e., N^}^ 
x I, where there are Nhii HII regions of size I. 

c 03-09. 

d Assuming a Salpeter IMF, m upP e r —80 Mq and m; o „, er = 1.0 

Mq. 
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Table 4 

The low- resolution model results for He 2-10 



Parameter 


Models I&II 


Model III 


Model IV 


n e (cm" 3 ) a 


500-1000, 


5000-50000 


2500-25000 




5000- 10 4 






size (pc) a ' 6 


30-40, 


2-10 


2-10 


EM (pc cm- 6 ) a 


3-8 
2xl0 7 


5xl0 7 -2xl0 9 


~10 8 




>5xl0 8 






Mm i (M s ) 


3xl0 5 ,~10 4 


7000-40000 


10 4 



Solutions for minimum N i!/C 



N™y™ (s; 1 ) 


2.5xl0 52 


4xl0 52 


5xl0 52 


n e (cm 3 ) 


5000- 10 4 


5000-50000 


5000-10 4 


size (pc) b 


4-7 


5-10 


3-7 


EM (pc cm" 6 ) 


>5xl0 8 


~10 8 


3xl0 8 


o /Qcont 
Oth/ & obs 


0.2-0.5 


0.3-0.7 


0.3-0.6 


Mhii (Mq) 


-10 4 


2xl0 4 


10 4 -50000 


No. O stars c 


1825 


2900 


3650 


Mass of stars (M ) d 


3.5xl0 5 


5.5 xlO 5 


7xl0 5 



a The two rows correspond to two classes of models, both of which 
are consistent with existing observations. The models predict dif- 
ferent RRL and continuum flux densities at other frequencies (see 
Figure 3) 

b 'Effective' sizes are quoted for Models III and IV, i.e., x I, 

where there are Nmi HII regions of size I. 

c 03-09. 

d Assuming a Salpeter IMF, m upper =80 M Q and m; ou)er =1.0 M©. 



Table 5 

Model results for 8.3 GHz line and the high resolution continuum data of KJ99 and TBHOO 



Parameter NGC 5253 He 2-10 



n e (cm 3 ) 


5xl0 4 -10 6 


5000-10000 


size (pc) 


1.0-2 


-7 


Nl, c (s- 1 ) 


>(2.5-4.0)xl0 52a 


>2.5xl0 52 


T c 


>10 


-0.6 


[SH92a]model / [Sff92a] obs 


0.1-0.35 


-1 


M HII (M Q ) 


>2000 


>2xl0 4 


No. stars b 


> 1900-2900 


>1900 


Mass of stars (M ) c 


>(3.5-5.5)xl0 5 


>3.5xl0 5 



a The range 2.5-4 is over various models. 
b 03-09. 

c Assuming a Salpeter IMF with m upper =80 M and 
mi ower =1.0 M Q . 



